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Fundamental principles of mapping 3-dimensional quasi- 
particle dispersions in the valence band using angle-resolved 
photoemission spectroscopy are discussed. Such mapping is 
intrinsically limited in accuracy owing to damping of the fi- 
nal states, resulting in equivalent broadening in the surface- 
perpendicular wavevector. Mechanisms of the intrinsic accu- 
racy are discussed in depth based on a physically transparent 
picture involving interplay of the final- and initial-state spec- 
tral functions, and illustrated by photoemission simulations 
and experimental examples. Other interesting effects of 3- 
dimensional dispersions include 'ghost' photoemission peaks 
outside the Fermi surface and finite peak width at the Fermi 
level. Finally, optimization of the experiment on the intrinsic 
accuracy is discussed. 

Keywords: Bandstructure; Photoemission; 

3-dimcnsional effects; Final states; Intrinsic effects 

I. INTRODUCTION: WHY INTRINSIC 
ACCURACY? 

The quasiparticle bandstructure E(\t), which reflects 
the peaks of the spectral function A(E, k) depending on 
energy and wavevector k, is the key property of the crys- 
talline solids. Angle-resolved photoemission (PE) spec- 
troscopy (for reviews see, e.g., [1,2]) is the main experi- 
mental method to map E(k) with resolution in energy E 
and, in principle, in the 3-dimensional k {3D band map- 
ping) . This is based on conservation of k in the process 
of photoexcitation from the occupied initial state to un- 
occupied final state in the bulk of the crystal. However, 
the PE experiment is inherently performed on the crys- 
tal terminated by the surface. By virtue of the remaining 
surface-parallel periodicity, the surface-parallel wavevec- 
tor k|| is conserved in the whole system and can be deter- 
mined from the surface-parallel vacuum wavevector K|| 
with the accuracy limited only by instrumental resolu- 
tion. As the surface-perpendicular periodicity is however 
broken, the surface-perpendicular wavevector k± is not 
conserved in the whole system and can not be directly 
determined. Control over k± is therefore a fundamental 
problem of 3D band mapping. 

The control over k± has two aspects. First, k± is dis- 
torted when the photoelectrons exit from the bulk of the 
solid into vacuum. As k± is conserved in the very pho- 
toexcitation process in the bulk, it can nevertheless be 
recovered if the fc^-dispersion of the final state E$ (k±) or 
that of the initial state E' l (k±) back in the bulk is known. 



Commonly an empirical frcc-clcctron-like approximation 
is used here for the final states, though the accuracy 
of this can often be insufficient. Nevertheless, the true 
quasiparticle E' (k±) including the non- free-electron and 
self-energy effects in the final state can be determined us- 
ing very-low-energy electron diffraction (VLEED) spec- 
troscopy (see [3-5] and references therein) based on that 
the PE final states, as proved by the one-step theory of 
PE, are the time-reversed LEED states [6,7]. It is a com- 
mon point of view that knowledge of the final-state E(k) 
exhaust the whole problem of 3D band mapping. Indeed, 
variation of the final-state energy Ef in the PE exper- 
iment (commonly achieved by variation of the photon 
energy hv) would through known E^(k±) be translated 
into variations of k±, and accurate 3D mapping of the 
occupied -E'(k) would be performed. 

There is another aspect however. Due to inelastic ab- 
sorption and elastic reflection from the crystal poten- 
tial the PE final state is intrinsically damped towards 
interior of the solid [6,7]. Such a confinement in the 
surface-perpendicular coordinate is equivalent, by the un- 
certainty principle, to certain intrinsic broadening in k±. 
The PE signal via this final state is then formed by all 
initial states within the k± -broadening interval, and the 
PE peak reflects an average of the quasiparticle E z (k±) 
in the valence band. A difference of the measured aver- 
aged /cj_-dispersion from the true E l (kj_) limits intrinsic 
accuracy of 3D band mapping. It is borne by fundamen- 
tal physics of the PE process, and can not be improved 
instrumentally. 

The intrinsic accuracy appears thus as an extrinsic fac- 
tor to the excited-state self-energy effects as deviations 
of the quasiparticle excited-state E(k), measured in the 
ideal PE experiment with undamped final states, from 
the ground-state bandstructure. 

The intrinsic accuracy problem has in pieces been tack- 
led in quite a few papers. Here, my aim is to delineate 
an entire concept, including essential physics and mecha- 
nisms of the intrinsic accuracy, and unveil its non-trivial 
effects encountered in the PE experiment. The concept 
is built upon a simple picture involving interplay of the 
final- and initial-state spectral functions, which gives a 
clearer physical insight compared to the explicit one-step 
PE theory with its heavy computational machinery. Fi- 
nally, I discuss optimization of the PE experiment with 
respect to the intrinsic accuracy of 3D band mapping. 
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II. BASIC PHYSICS OF THE PE PROCESS 

In the following analysis of the PE process it is implied 
that E(k) is the quasiparticle bandstructure, already in- 
corporating the excited-state self-energy effects. We will 
neglect the surface effects and concentrate on the bulk 
derived states with significant dispersion in k± (so-called 
3D states). The ideal E- and k|| -resolution of the PE 
experiment is assumed. 



A. Initial and final states of PE 

First we will recall physics of the PE initial and final 
states based on the ideas from [1,2,8]. For simplicity we 
will assume that their wavefunctions have only one Bloch 
wave constituent. 

The final-state wavefunction is the time-reversed 
LEED wavefunction. As sketched in Fig.l (the upper 
panels represent the final state), the wavefunction is de- 
scribed by a Bloch wave which is damped in the surface- 
perpendicular direction rj_ towards the crystal interior 
[3,5,9,10]. There are two sources of such damping: 1) 
inelastic absorption, expressed by absorption potential 
Vi connected with the electron lifetime as 2Vi = h/r e , 
and 2) elastic reflection from the crystal potential in the 
final-state band gaps. The damping is described by com- 
plex k^j_ component of the final-state wavevector k-^, with 
Imk ± connected to the photoelectron escape length A as 
2Imk^ = A -1 (the factor of 2 comes from squaring the 
wavefunction amplitude for intensity). 

The corresponding final-state dispersion E^(k±), 
shown in Fig.l by black lines, is radically different from 
this one if the Bloch waves were undamped (propagat- 
ing), as shown in gray: the dispersions as a function 
of Rek± are smooth and gapless, passing continuously 
through the band gaps of the undamped E(k) [3,5,9,10]. 
Such regions will in the following be somewhat loosely 
referred to as final-state band gaps. Note that due to 
surface-parallel invariance of the PE process the Bloch 
waves are undamped in the surface-parallel direction, re- 
taining real k|| and gapped k|| -dispersions [5]. 

The final-state spectral function A^(E,k), due to the 
finite photoelectron lifetime, is in principle characterized 
by certain distribution in E. However, in the PE experi- 
ment this is reduced to the 5-function at the energy given 
by the PE analyser. Due to damping of the final-state 
wavefunction A?(E, k) is characterized, on the contrary, 
by Lorentzian distribution in real k±, as shown in Fig.l, 
centered at kj_ = Rek^ of the final-state Bloch wave and 
having the fullwidth 5k± = 2Imk^ (in the angle-resolved 
PE experiment k|| is fixed and can be omitted). There- 
fore, Af(E, k) is given by 
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where the nominator stands for the normalization 
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J Al (kj_)dk^ = 1. The final state of PE is thus char- 
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acterized by fixed E and broadening in k± whose width 
is determined by Imk^_ of the damped final-state Bloch 
wave [1,2,8]. 

Nature of the PE initial state is different and com- 
plementary. The initial-state wavefunction is basically 
undamped: Removal of holes by the finite hole lifetime 
is compensated by their simultaneous generation by the 
electromagnetic field, which is almost homogeneous, ex- 
tending into the crystal interior over a light absorption 
length of ^100 A, very large compared to the final-state 
damping length. As sketched in Fig.l (the lower panels 
represent the initial state), initial-state wavefunction is 
therefore described by an almost undamped Bloch wave, 
having well-defined real k± . 

The corresponding initial-state dispersion E t (k±) 
shows up sharp and gapped /c^-dispersions typical of the 
undamped states, radically different from the final-state. 

The initial-state spectral function A t (E,'k) is charac- 
terized, owing to well-defined k±, by distribution in real 
k± reduced to the <5-function. Due to finite hole lifetime 
t/j it is characterized, on the contrary, by Lorentzian dis- 
tribution in energy, as shown in Fig.l, centered at the 
band energy E l (k±) and having the fullwidth SE = h/rh- 
Therefore, A l (E,k) is given by 
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where the nominator ensures J A l (E)dE = 1. The ini- 
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tial state of PE is thus characterized, complementarily to 
the final state, by fixed k± and broadening in E whose 
width is determined by the hole lifetime [1,2,8]. 

Note that some implementations of the one-step PE 
theory [7,11] describe the hole lifetime effects by damp- 
ing of the initial states. However, such an approach is in 
clear contradiction with the PE experiment, which unam- 
biguously finds in the valence band the gapped E l (k±) 
dispersions typical of undamped states. The same funda- 
mental flaw is implied by an PE lineshape analysis (see, 
e.g., [12,13]) where the i?-broadening of the initial state is 

replaced by an equivalent k± -broadening ~ SE (^§§j^ 
Although such a replacement does not alter the lineshape 
while the k± -dispersions remain linear, it leaves out some 
of the intrinsic effects discussed here. 

Inverse photoemission (IPE) is commonly treated as 
the time-reversed PE process [14] in which the IPE final 
state in the conduction band is equivalent to the PE ini- 
tial state in the valence band. However, the above argu- 
ments suggest a fundamental difference between the two: 
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as the IPE final state is excited by the damped LEED 
wavefunction, it is, in contrast to the PE initial state, also 
damped. Correspondingly, the conduction band E(k) 
measured in the IPE experiment should be characterized 
by gapless k± -dispersions similar to the PE final-state 
dispersions in Fig.l. Surprisingly, such a fundamental 
difference between the IPE and PE processes, tracing 
back to fundamental principles of quantum mechanics, 
has so far been addressed nor theoretically neither ex- 
perimentally 

B. Development of the PE spectrum: Intrinsic 
fcx-resolution 

We will now apply the above principles to analyze the 
PE response of 3D states. We will use a simple formalism 
based on interplay of the spectral functions [1,5,8], which 
is extremely simple and gives a clear physical insight. Its 
accuracy, despite principal limitations compared to the 
explicit one-step PE theory (see HIE), is nevertheless 
sufficient for evaluation of the PE peak profiles [15]. 

As illustrated in Fig. 2, the photoemission current 
I(Ef,E l ) at the final-state energy E$ and initial-state 
energy E % — E^ — hv is formed by adding up the elemen- 
tary photocurrents dI{E^ \E l ) from all dk± intervals, in 
which direct transitions take place, with k± spanning the 
whole Brillouin zone (BZ). Each dI(E^ , E l ) is weighted 
by the above spectral functions, A^(k±) to express close- 
ness of k± to kj_ of the final-state Bloch wave, and A l (E l ) 
to express closeness of E' 1 to the initial-state band en- 
ergy E l (k±) at this k± with SE varying along E l (k±). 
Moreover, each dI{E^E % ) is multiplied by amplitude 
factors of the final-state surface transmission |T^| 2 and 
photoexcitation matrix element |M/i| . This yields 
dI{Ef,E l ) = dk ± \Tt'\ 2 \M fl (k ± )\ 2 Af (k ± )A l (E l ). To 
obtain the whole photocurrent, this expression should be 
integrated in k±: 

HE^E^cx J dk ± \Tf\ 2 \M fl (k ± )\ 2 ■ (3) 

— OO 

5k± 

(fc ± -fc«) 2 + (<5fc ± /2) 2 ' 
SE 

(E* - E i (k±)f + (SE/2) 2 

This integral can be evaluated at negligible computa- 
tional cost (it can even be performed analytically using 
Taylor expansion of E l (k±) [5]). 

Cuts of I(E^,E l ) along the E f — const line cor- 
respond to the constant-final-state (CFS) measurement 
mode, along the hv = E^ —E l = const line to the energy- 
distribution-curve (EDC) mode, and along E % — const to 
the constant-initial-state (CIS) mode. Fig. 2 illustrates, 



for simplicity, the PE spectrum taken using the CFS 
mode in which k\ is constant throughout the spectrum. 

The PE spectrum shows up, whichever the measure- 
ment mode, a peak centered roughly at the energy 
E l (kP L ) dictated by the direct transition at k\ between 
the final and initial bands. Broadening of the peak re- 
sults from the final-state fcj_-broadening combined with 
the initial-state _E-broadening. In this context the k±- 
broadening appears as an intrinsic k ^-resolution of the 
PE experiment (in the following we will use these terms 
interchangeably depending on the context). Note that 
this resolution, in contrast to the instrumentally limited 
resolution in kii, is limited by intrinsic mechanisms in- 
volved in the PE process and can not be improved in- 
strumentally. 

C. Regimes of the PE experiment 

If the fcj_-resolution is at its best extreme 5k± — ► 0, 
then Af (k±) becomes (5-function d(k± — k^_) and the in- 
tegral (3) for I{E* , E l ) is reduced to A 1 (E) whose max- 
imum points to the initial band energy E l (kj_) in the 
k-point dictated by the direct transition. The PE peaks 
will then exactly follow the 3D quasiparticle dispersions 
in the valence band. In practice, such a situation is real- 
ized under the condition 

5k ± < kf z (4) 

where k^ z is the surface-perpendicular dimension of the 
BZ. This condition identifies the so-called bandstructure 
regime (BS-rcgimc) of the PE experiment which ensures 
accurate 3D band mapping [2,6]. In the real space, it 
reads as A 3> c_i_, where c± is the surface-perpendicular 
dimension of the unit cell. 

If the PE final state comprise a few Bloch waves, the 
condition (4) should be generalized as 5k ± <C where 
N is the number of the corresponding final bands. Of- 
ten in the literature (see, e.g., [6]) this number is con- 
fused with the number of all unoccupied bands available 
for given E^ and k|| ; as the latter becomes with energy 
progressively immense, an erroneous conjecture is made 
that at high energies any 3D band mapping is impos- 
sible. However, in the multitude of the available unoc- 
cupied bands only those are effective as the final bands, 
whose Bloch waves provide effective transport of the pho- 
toelectrons out of the solid by significant coupling to the 
outgoing plane wave [5,16]. Typically there are only a 
few (N = 1 — 2) final bands in this sense (see examples 
in III A). 

If the A:_i_-resolution is at the opposite extreme Sk± — > 
oo, any resolution in k± is lost, A? (k±) becomes con- 
stant and the integral (3) yields, assuming that |T^| 2 
and |M/i| 2 are independent of k± and E t (k±) is linear, 
the 1-dimensional DOS (1DOS) f|t. The PE peaks will 
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then manifest the 1D0S maxima independently of k^_. In 
practice, such a situation is achieved under the condition 
Sk± > k^ 2 , which identifies the so-called IDOS-regime 
[2,6]. 

Interpretation of the PE data between the bandstruc- 
ture and 1DOS regimes is doubtful [6,17]: the PE peaks 
retain some /i^-dispersion, although damped by strong 
averaging in k± , but they can be far from the true quasi- 
particle E % (k±). 

Evolution of the regimes of the PE experiment with 
the fcj_-resolution will now be illustrated by a PE sim- 
ulation. The model employs the initial state whose pa- 
rameters are chosen to resemble roughly the sp-band of 
Cu [5]: It is described by the free-electron dispersion 
E l (k±) = ^—k\ + Vboo an d the energy broadening 8E 
varying linearly as a function of energy to match zero 
at the Fermi level Ep and an experimental value of 1.7 
cV at the bottom of the band. The Fermi-Dirac cut is 
introduced as the ideal step function. The photocur- 
rcnt I(Et ,E l ) is evaluated by integration of the spec- 
tral functions in k± (3), assuming constant amplitude 
factors |T^| 2 and |ikf^^ | 2 . The integration also extends 
above the Fermi vector k^_ to include certain contribu- 
tion from the tails of A l {E) in the unoccupied region of 
the valence band (see IIIB). For clarity, the PE spectra 
were evaluated for the CFS mode, in which k\_ is con- 
stant (deviations from the EDC mode are insignificant, 
because normally the final-state is large enough to 
ensure small variations of k\ through the peak profile). 

The simulation was performed for four Sk±^ values: zero 
for the ideal fcj_-resolution, 0.2 \TX\ as a representative 
for the BS-regime in a real PE experiment, 0.6|rX| for 
degraded fcj_-resolution, and 2|rX| for the IDOS-regime. 
The results are shown in Fig. 3 as series of the PE spectra 
with k\ scanning along TX (upper panels, each series 
in different intensity scale) and positions of the spectral 
peaks mapped on top of the true valence band (lower 
panels). 

In the ideal k± -resolution limit the PE peak reflects the 
initial-state A % (E) and ideally follows the true E % (k\ ). 
Interestingly, upon passing fc^ there appears an ex- 
tremely small peak just below Ep\ it replicates the low- 
energy tail of the spectral function A 1 (E) whose maxi- 
mum is already above Ep (see IIIB). 

In the BS-rcgimc (5k± = 0.2|rX| in our simulation) 
the /^-resolution remains sufficient for the PE peaks to 
closely follow E' l (k±). This regime can be used for ac- 
curate 3D band mapping. Two peculiarities should be 
noted here: 1) Despite A l (E) becomes the 8- function 
upon approaching Ep, the PE peak does not sharpen 
up to singularity; 2) When kj_ passes k± and enters into 
the unoccupied region, the peaks do not disappear but 
become 'ghost' peaks which are highly asymmetric and 
have their maximum stationary in energy just below Ep. 
Both peculiarities are appearances of the k± -broadening 



and will be analyzed in IIIB and IIIC. 

As the k± -resolution degrades (Sk± = 0.6|rX|), the 
concomitant averaging of E l (k±_) causes asymmetry of 
the peaks and strong intrinsic shifts of their maxima from 
the true E l {k±). Although the PE peaks deceptively 
keep dispersion on k° ± , their use for 3D band mapping 
will return a distorted valence band dispersion. Note ap- 
pearance of an additional shoulder-like structure, which 
originates from the large 1DOS in the bottom of the band 
and becomes notable as the main peak shifts upwards. 

In the IDOS-regime (Sk± = 2|rX|) the fc^-resolution 
is completely lost. The maxima of the PE peaks be- 
come stationary in energy reflecting the maximum of the 
valence band 1DOS (and a tiny replica of the A l {E) sin- 
gularity at Ep, omitted from the plot). The kj_ changes 
show up only in the intensity modulation. Note how- 
ever that the PE peaks are notably shifted from the true 
1DOS maximum. 

The regimes of the PE experiment in general follow the 
well-known "universal curve" X(E) [1,2,6] which gives the 
k± -resolution as 5k± = A -1 . The BS-regime holds while 
A remains large, which takes place at low E? below or 
not far above a plasmon excitation energy of 20-30 eV 
(see examples in III A) . Upon further increase of energy 
A reaches its minimum around 50-100 eV, and the PE ex- 
periment enters into the IDOS-regime, although in prac- 
tice the fcj_-dispersion signatures are never completely 
suppressed. The BS-regime recovers in the high-energy 
region above ~ 300 eV where A rises again (despite at 
these energies the unoccupied bands become immense in 
number, only one of them remains effective in the PE fi- 
nal state in the sense of effective coupling to the outgoing 
plane wave [16]). 

The high-energy region is characterized by low PE in- 
tensity resulting from small valence band crossection and 
from the Debye- Waller factor, in which the increase of hv 
reduces the PE intensity on equal footing with temper- 
ature. However, the advent of new high-brilliance syn- 
chrotron radiation sources and multidetection PE analy- 
sers has allowed for PE experiments in this region with 
reasonable count rates and resolution in E and k|| [18,19]. 
k-resolved experiments were reported, for example, in 
[20] where the sp-band of Al was found to be mirrored in 
the ^-dispersion of the PE signal up to 750 eV. Inter- 
estingly, this study has found an increase of the surface 
to bulk signal intensity ratio with hv, opposite to the 
tendency expected from the " universal curve" , and at- 
tributed this surprising phenomenon to intrinsic phonon 
excitation effects involved in the PE process. These re- 
sults call for further studies on the role of the phonon 
excitation in high-energy PE spectroscopy. 
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III. MECHANISMS LIMITING THE INTRINSIC 
ACCURACY 

Degradation of the fc^-resolution, as we observed 
above, results in intrinsic shifts of PE peaks from the 
true valence band /c^-dispersions, limiting the intrinsic 
accuracy of 3D band mapping. These shifts can originate 
through various mechanisms, which will now be analyzed. 

A. Non-linearity mechanism 

Degradation of the fc^-rcsolution can crawl into the 
position of the PE peaks, firstly, by non-linearity of the 
valence band fcj_-dispersion. As evident from Fig. 2, in 
this case the integral number of the dk± -intervals within 
the k± -broadening will be different for the initial states 
whose E l is above the direct-transition energy £^(fc°) 
and for the states whose E' 1 is below E l (kP L ), being larger 
where the 1DOS is larger. Correspondingly, the PE peak 
becomes asymmetric. Its maximum somewhat shifts to- 
wards larger integral number of dk±, deviating from the 
true E l (k^). In the band interior the PE peaks shift to- 
wards the band edges, where the 1DOS is larger. Near 
the band edges the peaks, on the opposite, shift towards 
the band interior (so-called in-band shifting [4,5]) be- 
cause beyond the band edge there are no states and the 
1DOS is zero. We will refer to this mechanism to limit 
the intrinsic accuracy as the non-linearity mechanism. It 
was behind the intrinsic shifts observed in the PE simu- 
lation in Fig. 3, where the PE peaks clearly followed the 
characteristic pattern of shifting towards larger 1DOS. 
Note that the intrinsic shifts due to the non-linearity 
mechanism depend, apart from 5k± and non-linearity of 
E l (k±), on the initial-state SE. 

The non-linearity mechanism can be illustrated by ex- 
perimental data for VSe2 (for a detailed discussion see 
[4]). Information about the final-state dispersions and 
lifetimes, required for determination of Sk± , was achieved 
in this case by VLEED (see IV). The experimental en- 
ergy dependence of Vi is shown in Fig. 4 (a). It is char- 
acterized by a sharp increases of Vi at the plasmon ex- 
citation energy hu! p . The experimental E^(k±) is shown 
in Fig. 4 (b). Note that these final band can hardly 
be described by free-electron-like dispersions: there are 
two bands, and each can be free-electron fitted only lo- 
cally with the inner potential strongly depending on E> 
and k|| (such a complicated structure of the final states 
is in fact typical of many materials with large unit cell, 
containing more than one atom, due to backfolding of 
many free-electron bands into the first BZ and their mu- 
tual hybridization). An estimate of bk\_ was obtained, 
neglecting some increase in the final-state band gaps, 

as 2Vi (ffj-) (see IV). The obtained Sk±_ values are 
shown in Fig. 4 (6) superimposed on the final bands. 



Following the increase of V{, the fc^-resolution degrades 
above hu) p . The results of PE band mapping, the binding 
energy of the PE peaks as a function of fc", are shown 
in Fig. 4 (c) with distinction for the final states below 
and above thu p . Near the band edges of the 3D bands the 
experimental points above Huj p are systematically shifted 
towards the band interior compared to those below hw p . 
It should be stressed that here we have an explicit exam- 
ple how the same initial states yield different energies of 
the PE peaks depending on the fc^-resolution. The ob- 
served shifts demonstrate the intrinsic effect of in-band 
shifting, occurring mainly by the non-linearity mecha- 
nism. Interestingly, increase of towards 50 eV results 
in further increase of Sk± and suppression of any disper- 
sion of the PE peaks [21]. 

Our analysis demonstrates that for VSe2 the energy 
region of accurate 3D band mapping is exhausted by low 
E* below or not far above Huj p . Such a situation is in fact 
fairly general, because a pronounced plasmon threshold 
is typical of many materials with some exceptions such 
as Cu and Ni. The intrinsic accuracy effects similar to 
VSe 2 were also observed for TiS 2 [4], WSe 2 [22], Ge [17], 
etc. 

In the final-state band gaps the Bloch waves experi- 
ence additional damping due to elastic scattering from 
the crystal potential [5,23]. Although this is normally 
much weaker than the damping due to Vi, as in the 
above example, for materials with exceptionally large 
final-state band gaps the additional damping can cause 
notable degradation of the fc^-resolution. Graphite is 
a typical example of this (for a detailed discussion see 
[24,25]). Its final bands were simulated using the empir- 
ical pscudopotential method based on the VLEED data 
for their energies in the T-point (the simulation is there- 
fore less accurate in the BZ interior) and Vi. The re- 
sulting Ef(k±) is shown in Fig. 5 (a) and (b) as a func- 
tion of Imk ± and Rek ± , respectively, using the doublc- 
BZ representation to reflect the dipole selection rules in 
graphite. The final states are characterized by a strong 
increase of Sk± = 2Imk^_ in the two band gaps in the 
T-point. Based on these data, we simulated the PE re- 
sponse of the valence 7r-band, having the 3D character, 
using the same formalism (3). The results are shown in 
Fig. 5 (b). When the final-state moves towards the T- 
point and enters the band gaps, progressive degradation 
of the fc^-rcsolution results, by the non-linearity mech- 
anism, in increasing of the in-band shifts. In the lower 
band gap, which has the maximal width, the increase of 
5k± is particularly strong; the concomitant in-band shift- 
ing is so dramatic that it overcomes the trend dictated by 
the band dispersions themselves and reverses the disper- 
sion of the PE peaks compared to the true 7r-band. Such 
an unusual reversed dispersion has indeed been found in 
the PE experiment [25]. 

Typically, intrinsic shifts near the band extrema due 
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to the non-linearity mechanism are of the order of a 
few tenths of eV. The experiment often suggests how- 
ever larger shifts [4,21,22]. As degradation of the k±- 
resolution goes along with reduction of the escape length 
A, this reveals certain contribution of surface effects such 
as suppression of the 1DOS in the band extrema (see 
HID). It should be noted that the condition (4) identi- 
fying the BS-regime is particularly restrictive for layered 
materials owing to their small k^_ z ( see [26] for a spe- 
cialized survey). 

B. 'Ghost' peaks at Ef 

An interesting phenomenon occurs when a 3D band 
crosses the Fermi surface. This situation is described by 
the PE simulation in Fig. 3: When the direct-transition 
k\ moves outside the Fermi surface and starts sampling 
the unoccupied part of the valence band, the PE peaks, 
mysteriously, do not disappear. Instead, they become 
asymmetric and stationary in energy just below the Fermi 
level Ep, and only gradually reduce in amplitude. Ex- 
perimentally, such a phenomenon has been observed, for 
example, in a famous PE study on Na, where the sp-band 
crosses Ep [27,28]. We will refer to this phenomenon as 
'ghost' peaks. It has in fact two contributions due to 
different mechanisms. 

The first contribution is due to a general mechanism 
taking place for 3D as well as 2D bands. It comes from 
the unoccupied states next to Ep: The spectral function 
A l (E), centered at the band energy above Ep, protrudes 
its low-energy tail, however small the remaining ampli- 
tude, below Ep and gives therefore some PE intensity 
(such an effect for 2D bands is discussed, for example, in 
[29,30]). Such a contribution can be isolated in our PE 
simulation by restricting the integration (3) to k± > k±. 
It survives even in the limit of the ideal fc_L-resolution, 
see Fig. 3. In any case, the 'ghost' peak intensity due this 
mechanism is very small. 

The second contribution, much predominating for the 
3D bands, comes from the occupied states: As illustrated 
in Fig. 6, even if fc° moves outside the Fermi surface, 
the dk±-sta.tes at the initial-state E % below Ep give cer- 
tain 'ghost' intensity by virtue of being accessible via the 
k± -broadening. Upon moving fc^ further from k^_ the 
(k±) profile shifts, reducing its amplitude at the ini- 
tial band and, concomitanly, the 'ghost' peak amplitude. 
Note that this mechanism gives a significant 'ghost' peak 
intensity even if the initial-state ^-broadening is vanish- 
ing which takes place upon approaching Ep. 

The 'ghost' peak phenomenon complicates identifica- 
tion of the Ep crossings in mapping of the Fermi surface 
formed by 3D bands. A well-known example is NbSe2, 
whose Ap z band forms a 3D pocket of the Fermi sur- 
face with extension in k± small compared to the k±- 
broadening; direct mapping of this pocket is precluded 



by non-dispersive 'ghost' intensity [31,32]. Even if the 
mapping is performed as a function of ky, the 'ghost' 
peaks stay at Ep through the entire interval of ky where 
at least some part of the 3D band remains occupied. It 
should be noted that often the PE data, especially for 
quasi-2D materials, is interpreted neglecting the 3D ef- 
fects; in this case such 'ghost' peaks in the ky -dispersion 
are misinterpreted as a sign of unusual self-energy effects. 
A way to circumvent the 'ghost' peak problem in the 
Fermi surface mapping was suggested in [5]: one varies 
ky , simultaneously changing E* to remain in the E' 1 (fcjj 
minimum. 

C. Linewidths near Ep 

Similar in origin to the 'ghost' peaks is an interesting 
PE linewidth behavior upon approaching Ep: As seen 
in our simulations in Fig. 3 for the BS-rcgimc, the PE 
peak from the 3D band, mysteriously, retains a finite 
linewidth up to Ep without sharpening up as might be 
expected from the behavior of the initial-state A l (E), 
which becomes singular at Ep because of the vanish- 
ing .E-broadening. The origin of such a counter-intuitive 
linewidth behavior is again due to the 3D character of the 
band and the k± -broadening [13,33]: As one can conjec- 
ture from Fig. 6, even if the i?-broadening is zero, the 
initial states around the direct-transition energy E t (k <: ^) 
are accessible via the k± -broadening and contribute to 
the PE peak width. In the ideal fc^-rcsolution limit, as 
demonstrates our simulation in Fig. 3, the PE peak well 
becomes singular towards Ep. 

An experimental example of such an effect of the k±- 
broadening on the linewidths can be seen, for exam- 
ple, in the PE peak from the sp-band of Cu whose 
linewidth increases upon approaching Ep [15]. In [13,33] 
it was shown that the same effect can explain an un- 
usual linewidth energy dependence near Ep observed in 
Bi 2 Sr 2 CaCu208, which otherwise would indicate devia- 
tions from the Fermi liquid picture. 

Significant PE intensity next to Ep, seen in all regimes 
of the simulation in Fig. 3, also originates from existence 
of a 3D band crossing Ep and the fcj_-broadening. Such 
a peculiarity is typical of the PE experiments on 3D ma- 
terials such as Cu. 

D. Other mechanisms 

Intrinsic shifts in the 3D bands can also appear due 
to a matrix element mechanism: If Mfi as a function of 
energy and k± undergoes sharp variations within the k±- 
broadening, a contribution of the states above and below 
E l (kP L ) is different (see Fig. 2). The peak becomes asym- 
metric and it maximum experiences an intrinsic shift 
from the true band energy. Such a mechanism has been 
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identified, for example, in TiS2 [4]. Similar effects due to 
should be less pronounced, because the energy vari- 
ations of are normally smooth and become sharper 
only at extremely low , as can be seen in VLEED 
data (see, e.g., [4,24,34]) 

The above mechanisms were all based on the bulk 
bandstructure picture, in which the bandstructure near 
the surface is considered identical to the bulk one. This 
implies that the crystal potential, and thus the wavefunc- 
tions, ideally repeat the bulk ones up to the surface. In 
fact, within a few atomic layers of the PE escape depth 
they can experience notable modification. Such surface 
effects also give rise to intrinsic shifts. For example, the 
1DOS singularities in the band extrema, characteristic 
of the bulk E(k), are smeared in the local DOS near the 
surface (see, e.g., [35-37]), which results in an increase 
of the in-band shifting. The surface effects are harder 
to control compared to the fc^-rcsolution effects. Qual- 
itatively, their contribution reduces with increase of A 
simultaneously with reduction of the intrinsic shifts of 
the bulk origin due to the concomitant improvement of 
the fc_L-resolution. 

Surface photoelectric effect, generated by abrupt 
change in the dielectric response at the surface barrier 
and described by the V ' A part of M/,, can influence the 
PE signal from the 3D bands through interference of the 
surface and bulk PE components (see, e.g., [38-40]). Nor- 
mally this effect only causes some asymmetry of the PE 
peaks. However, near the surface and bulk plasmon exci- 
tation energies, especially for free-electron metals such as 
Na, the lincshapc changes can be drastic and return large 
intrinsic shifts [28,41]. In [40] it was demonstrated that 
the surface and bulk components can be experimentally 
separated by changing the light incidence angle. 

E. Intrinsic accuracy within the one-step PE theory 

The intrinsic effects are naturally embodied in the one- 
step PE theory (see, e.g., [6,7]) which finds the PE inten- 



sity as a matrix element 



Ap 



between the 



final- and initial-state wavefunctions <I>^ and $ l . Inter- 
ference between different Bloch wave constituents in the 
final and initial states, matrix element effects and sur- 
face effects, all ignored in the simplified formalism (3) 
used above, are naturally included here. However, a 
heavy computational machinery often obscures the phys- 
ical mechanisms. 

An example of the intrinsic accuracy analysis using the 
one-step PE theory can be found in [42] in application to 
TiSe2. Similarly to the above results for VSe2, for the 
final states above hu> p the -resolution was found to 
degrade, causing strong broadening and intrinsic shifts of 
the PE peaks. Surface effects due to smearing of the bulk 
1DOS singularities have also been identified. Another 



example is a calculation on Na [28], which described the 
'ghost' peaks at Ep and in-band shifting in the bottom 
of the valence band. 



IV. OPTIMIZATION OF THE PE EXPERIMENT 
USING VLEED 

3D band mapping remains relevant if the intrinsic 
shifts are small. Of the factors affecting these shifts, 
the non-linearity of the fc^-dispersion, initial-state E- 
broadening and, partly, variations of M/j are all due to 
inherent properties of the valence bands under study and 
remain beyond control by the experimcntator. The k±- 
resolution, on the contrary, is an exclusively final-state 
property, and by tuning E* can taken to values which 
ensure negligible intrinsic shifts. This is achieved in the 
BS-regime Sk± <C k^ z (see II C). The exact energy 
ranges of this regime are however much material depen- 
dent. Knowledge of the fc^-rcsolution energy dependence 
is therefore required. 

VLEED, covering the energy range below ~40 eV, has 
recently been established as the experimental method 
giving the most direct access to the PE final states with 
resolution in the 3D wavevector (see [3-5] and references 
therein). In the VLEED spectra of elastic electron re- 
flection, the energies of the spectral structures give the 
critical points in the final bands. E* (k±), obtained from 
this experimental data, can further be used to control k± 
in the PE experiment. Broadening of the VLEED struc- 
tures, moreover, gives an estimate of the corresponding Vj 
values [4,34,43]. By performing calculations of the com- 
plex bandstructure with the experimental V [3,5,34], one 
obtains the fcj_-resolution as Sk± = llmk^ of the gener- 
ated Bloch waves (and simultaneously the energy depen- 
dence of A, as it was recently demonstrated for graphite 
in comparison with determination of A by the conven- 
tional ovcrlaycr method [43]). As such calculations can 
however be time-consuming, 5k^ can be estimated in a 
simplified manner directly from the V values and the 
final-state group velocity as [9,12,13] 



Sk ± ~ 2V, 



\dk±) 



(5) 



This estimate is accurate however only outside the final- 
state band gaps. Nevertheless, increase of 5k± associ- 
ated with additional damping in the gaps (see Fig. 5) is 
normally less significant, except for a few materials with 
exceptionally wide band gaps such as graphite. 

Use of the VLEED predictions on the fcj_-resolution 
allows the experimentator to optimize the PE experiment 
on the intrinsic accuracy by selecting the intervals 
which ensure the BS-regime (see the examples in III A). 

As a concluding remark, it should be noted that one 
can significantly improve accuracy of 3D band mapping 
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by computational modelling of the intrinsic accuracy ef- 
fects, as described above, and correcting the experimen- 
tal data correspondingly. This is often critical for proper 
evaluation of the true excited-state self-energy effects, for 
example in graphite [25] and Na [28]. 

V. CONCLUSION 

Fundamental principles of PE have been analysed, fo- 
cussing on emission from 3D states in the valence band. 
Damping of the final state due to inelastic and elastic 
scattering in the crystal results in equivalent broaden- 
ing in k±, which acts as the intrinsic fc^-rcsolution of 
the PE experiment. The PE peak is then formed as 
a matrix-element weighted average of the quasiparticlc 
fc^-dispersion in valence band, which can be shifted from 
the dispersion itself. Such shifts intrinsically limit the 
accuracy of 3D band mapping, which remains accurate 
only in the bandstructure regime, characterized by the 
intrinsic fc^-broadening small compared to the BZ ex- 
tension in k±. Upon degradation of the fcj_-resolution 
the PE experiment yields increasing intrinsic shifts and, 
finally, enters into the 1DOS regime. The intrinsic k±- 
resolution can be controlled using VLEED, which gives 
the final-state dispersions and lifetimes. The mechanisms 
forming the intrinsic shifts have been exposed, including 
non-linearity of the fc^-dispcrsion, variations of the ma- 
trix element and surface effects. Appearance of the in- 
trinsic shifts and transformation of the regimes of the PE 
experiment has been analysed using a physically trans- 
parent picture of interplay of the final- and initial-state 
spectral functions, and illustrated by PE simulations and 
experimental examples. The fc^-broadening gives rise to 
other surprising effects typical of 3D valence states such 
as 'ghost' peaks whose direct-transition k± falls outside 
the Fermi surface, and finite width of the PE peaks at 
the Fermi level. 
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FIG. 1. Character- 

istic behavior of wavefunctions, fcx-dispersions and spectral 
functions of the final and initial states of PE, shown, respec- 
tively, on the upper and lower panels. The final states are 
characterized by fixed E and broadening in k±, while the ini- 
tial state, complementarily, by fixed k± and broadening in 
E. 

FIG. 2. Development of the PE peak. Its broadening 
results from the final-state fcx-broadening, defining the in- 
trinsic fcx-resolution 8k±, combined with the initial-state 
iJ-broadening 8E. The maximum of the peak shows an intrin- 
sic shift towards the larger average 1DOS from the position 
E l (k A _) dictated by the direct transition at k±. 



FIG. 3. Simulated PE experiments on a model 3D band, re- 
sembling the sp-band of Cu, performed with different strength 
of damping in the final state. The corresponding 8k± val- 
ues of the intrinsic the fcx-resolution (indicated on top) are 
characteristic of the BS-regime (left), degraded fcx-resolution 
(center), and IDOS-regime (right). The results are shown as 
series of a few PE spectra (upper panels) whose offsets reflect 
scanning of the direct-transition fc° along FX, and the posi- 
tion of the spectral peak maxima on top of the true valence 
band E l (k±) (lower panels). Finite width of the PE peak 
upon approaching Ef and 'ghost' peaks remaining after the 
Ef crossing as both effects of the fcx-broadening. Degrada- 
tion of the fcx-resolution results in intrinsic shifts of the PE 
peaks, in this case is due to the non-linearity mechanism, from 
the true valence band. 

FIG. 4. Intrinsic accuracy for VSe2: (a) Energy depen- 
dence of the final-state Vt with a characteristic increase due 
to the plasmon excitation at hw v ; (b) Final-state bands (in- 
cluding, of all unoccupied bands, only those whose Bloch wave 
effectively couples to the outgoing photoelectron plane wave). 
The fcx-resolution is shown by shading; (c) Position of the 
PE peaks superimposed on DFT-LAPW calculations. The 
points obtained with E 1 below hui p are shown as solid dots, 
and above as open dots. A systematic in-band shifting in the 
3D bands, following degradation of the fcx-resolution above 
h~ujp, is mainly due to the non-linearity mechanism. 

FIG. 5. Intrinsic accuracy for graphite: (a,b) Empiri- 
cal pseudopotential simulation of the final bands based on 
VLEED data (black lines); the band gaps are emphasized 
by setting Vi to zero (gray); (c) Simulated position of the 
PE peaks (dense dots) from the valence 7r-band on top of 
its true dispersion. Severe degradation of the fcx-resolution 
Sk± = Imk ± in the final-state band gaps causes large in-band 
shifts due to the non-linearity mechanism, and even reverses 
the dispersion of the PE peaks in the lower gap. 

FIG. 6. Mechanism of the 'ghost' peaks from 3D bands 
crossing the Fermi surface. Even if the direct-transition fc° 
is in the unoccupied region, the PE signal comes from the 
occupied part of the valence band accessible via the final-state 
fcx-broadening. 
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